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The recent series of experiments on polarized lepton-nucleon scattering have pro- 
vided a strange new twist in the story of the nucleon, some of whose aspects are 
reviewed in these lectures. In the first lecture, we review some issues arising in the 
analysis of the data on polarized structure functions, focusing in particular on the 
importance and treatment of high-order QCD perturbation theory. In the second 
lecture some possible interpretations of the "EMC spin effect" are reviewed, princi- 
pally in the chiral soliton (Skyrmion) approach, but also interpretations related to 
the axial !7(1) anomaly. This lecture also discusses other indications from recent 
LEAR data for an ss component in the nucleon wave function, and discusses test 
of a model for this component. Finally, the third lecture reviews the implications 
of polarized structure functions measurements for experiments to search for cold 
dark matter particles, such as the lightest supersymmetric particle and the axion, 
after reviewing briefly the astrophysical and cosmological evidence for cold dark 
matter. 



1.1 Formalism 

The basis for our discussion will be the two spin-dependent structure functions 
Gi and G2: 



"Invited lectures at the Int. School of Nucleon Spin Structure, Erice, August 1995. 
''This section is updated from Ref. M. 



1 Polarized Structure EYinctions 



" dQ'^dv ^ 

Mn{E + E' cose)Gx{u, Q2) _ Q-^G2{v, Q^) 




In the parton model, these structure functions scale as follows in the Bjorken 
limit X = Q'^ /2Mnv fixed, oo: 

Mf,uGi{u, Q^) = gi{x, Q^) ^ gi{x) 

(1.2) 

MNiy''G2{u,Q^) = g2{x,Q^) ^ 52(2;) 

We will discuss the scaling structure function 92 later on, focusing for now 
on gi , which is related to the polarized quark distributions by 

9ii^) = IYI ^9[5T(a=) - nix) + q^{x) - qiix)] (1.3) 

for comparison, the unpolarizcd structure function F2 is given by 

F2{x) =Y^elx[q^{x) + qi{x) + q^{x) - qi{x)] (1.4) 

so that the polarization asymmetry Ai may be written as 

A, = ^iZi^^ (1.5) 
(^1/2 + ('■3/2 

in the Bjorken limit, where ai/2 and 0-3/2 are the virtual photon absorption 
cross sections. We will discuss later the dependences of the above formulae, 
as well as the transverse polarization asymmetry. 

Much of the interest in the polarized structure function gi is due to its 
relation to axial current matrix elements: 

(pl^^b) = {p\<iit^i5q\p) = {p\qRiiJ.qR - qLiij.qL\p) = 

(1.6) 

where qL,R = 1/2(1 T 75)?) Sfj, is the nucleon spin four-vector, and 

Aq= [ dx[q^{x)-q^{x) + q^{x)-qi{x)] (1.7) 
Jo 



Of particular interest is the matrix element of the singlet axial current 

E Wsg: {p\Al\p)= ^1-S,{P) (1.8) 

q—u,d,s q—u,d,s 

which is related in the parton model to the sum of the light quark contribu- 
tions to the proton spin. Prior to the series of measurements of polarized deep 
inelastic lepton nuclear scattering, information was available from charged cur- 
rent weak interactions on some axial current matrix elements. For example, 
neutron beta decay and strong isospin symmetry tell us thatl3 

Au - Ad = + = 1.2573 ± 0.0028 (1.9) 

and hyperon beta decays and flavour SU{3) symmetry tell us thati 

+ - ^ = ^ - 0.34 ± 0.02 (1.10) 

v3 v3 v3 

(for a recent discussion of the applicability of S U (3) symmetry see Ref. and 
references therein). Equations ( p_.g| ) and ( 1.10 ) give us two equations for the 
three unknowns Au, Ad and As. In principle, a third piece of information was 
availabled'EI in 1987 from nejutral current weak interactions. Measurements of 
vp and vp elastic scatteringlZl indicated that 

As = -0.15 ±0.09 (1.11) 

but this information was not generally appreciated before the advent of the 
EMC data discussed belowj^ At present there is a new neutrino experiment 
under way at Los Alamos u which is expected to significantly improve the 
precision of ( [l . 1 1[ ) (sec Rcf. |^ for a recent in-depth analysis). 

In the naive parton model, the integrals of the gi structure functions for 
the proton and neutron 



Jo 



(1.12) 



rUQ') ^ / dx g^x 







are related to combinations of the Aq. 



r;.l(lA„ + lAd + lA. 



(1.13) 
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The difference between tfie proton and neutron integrals yields the cele- 
brated Bjorken sum ruleliO 

r?(Q') - rr(g2) = ^(Au - Ad) x (i - a^mM + ■■■ (i-W) 

It is not possible to derive individual sum rules for F^'" without gjUpple- 
mentary assumptions. The assumption made by Ellis and Jaffe in 1973EJ was 
that As = 0, on the grounds that very possibly there were a negligible number 
of strange quarks in the nucleon wave function, and if there were, surely they 
would not be polarized. With this assumption, it was estimated that 

/ dxgP,{x, Q^)= — {AAu + Ad) (1 - ajir + ...) = 

(1.15) 

= 0.17±0.01 

It should be clear that this was never a rigorous prediction, and was only 
intended as a qualitative indication to experimentalists of what they might 
find when they started to do polarized electron proton scattering experiments. 

Rertpjbative QCD corrections to the above relations have been calcu- 
latedlia-y : 

[\gfix,Q')-gnx,Q')] 
Jo 

f{x) = l-x- 3.58x2 

and 

Jo 

= (^±^ \9a\ + ^ as) fix) + ^ AS(g2) h{x) : (1.17) 
h{x) = l-x - 1.096x2 - . . . 



9a\ fix) 



1 
6 

20.22x3 



(1.16) 



where x = a^iQ'^) ji:, and the dots represent uncalculated higher orders of 
perturbation theory, to which must be added higher-twist corrections which 
we will discuss later. The coefficients in ( |1.16D ,( I.17 ) are for Nf=3, as relevant 
for the Q2 range of current experiments. 
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With these corrections, the Bjorken sum rule is a fundamental prediction 
of QCD which can be used, for example, to estimate a value for as{Q^)- On 
the other hand, the individual proton and neutron integrals can be used to 
extract a value of As. 



1.2 The Helen of spin 



Early 
ments 

large errors. Over a 100Q_^theoretical and experimental papers were launched 
by the 1987 EMC result 



n polarized electron-proton scattering from SLAG- Yale experi- 
61111 "v^rere compatible with the prediction of equation (1.15) within 



gl{x, Q^) = 0.126 ± 0.010 (syst.) ± 0.015 (stat.) 
at (Q2) = 10.7 GeV^ 



(1.18) 



which was in prima facie disagreement with the dynamical assumption that 
As — 0. It is worth pointing out that the small- a; behaviour olg\{x) was crucial 
to this conclusion. The earlier SLAG- Yale data had large extrapolation errors, 
and the EMG data indicated behaviour different from that in simple dynamical 
models. They were, however, consistent eI with the naive Regge expectation t3 



(1.19) 



were the Q!i(0) are the intercepts of axial vector Regge trajectories which are 
expected to lie between and -0.5. A fit to the EMG data gaveB 



p 
9i 



-0.07I^;32 for x<0.2 



(1.20) 



for X < 0.2. 

Using equations ( |l.9| ) , ( I.IC ) , ( 1.18 ) and the leading order perturbative QGD 
corrections in equation (1.17) it was estimated 123 that 



Au 
Ad 
As 



0.78 ±0.06 
-0.47 ±0.06 
-0.19 ±0.06 



(1.21) 



Strikingly, these determinations corresponded to a total contribution of quarks 
to the proton spin 

AE = Au + Ad + As = 0.12±0.17 (1.22) 
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which was compatible with 0. 

This has sometimes been called the "pro ton spin crisis" , but we think this 
is an over-reaction. The result equation (1.22) was certainly a surprise for naive 
models of non- perturbative QCD, but it was not in conflict with perturbative 
|QCD. Moreover, shortly after the first data became available it was shown 
c3 that AE = occurs naturally in the Skyrme model, which is believed to 
reproduce the essential features of QCD in the large- A^c limit. Alternatively, it 
was suggested that the U{1) ptiaL anomaly and polarized glue might provide 
an alternative interpreta±iiDii£il or a significant suppression of the QCD 
topological susceptibilityEj c3 might play a key role, which would modify the 
nai've quark model predictions. These interpretations are discussed in more 
detail in Lecture II. 

In 1972 Richard Feynman wrote ". . . its [the Bjorken sum rule's] verifi- 
cation, or failure, would have a most decisive effect on the direction of future 
high-energy physics" . On the other hand, we think that the verification, or fail- 
ure, of equation ( 1.22 ) has only an indecisive effect, though a very interesting 
one. 



1.3 Evaluation of integrals 

Before discussing the interpretation of more recent data on polarized structure 
functions, we first review a few points that arise in the evaluation of the inte- 
grals r^'". It should not be forgotten that the QCD versions of the sum rules 
are formulated at fixed Q^. A generic deep-inelastic sum rule in QCD reads 



r(Q2) = Fo 



Ji>l 



as(Q') 



m>l ' 



(1.23) 



where Foo is the asymptotic value of the sum rule for 



the c„ are the 



coefficients of the perturbative corrections, and the dm are coefficients of the 
so-called mass and higher-twist corrections. On the other hand, the data are 
normally obtained at values of that increase on the average with xsj as 
seen in Fig. pTl] . 

It is therefore necessary for each individual experiment to interpolate and 
extrapolate to some fixed mean value of , as indicated by the dashed hor- 
izontal line in Fig. |l.l| . The quantity measured directly is the polarization 
asymmetry (1.5), which seems experimentally to have only small dependence 
on Q^. In particular, a recent analysis by the E143 collaboration EZl sees no 
significant dependence in Ai for 1 GeV^ , and the dependence seen 

at lower is compatible with the magnitude of the higher-twist correction 
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4 



(GeV^) 



<Q'> 



0.0 0.5 1.0 

Figure 1.1: In any given polarized lepton-nucleon scattering experiment, the range of 
probed is different in different bins of the Bjorken variable xbj ■ 



discussed later. Therefore experiments often assume that Ai is a function of 
X only, and then estimate 

^ A,ix,Q^)F2ix,Q^) A,ix)F2ix,Q^) 
^'^^'^ ^ 2x[l + i?(a;,Q2)] ' 2x[l + Rix , Q^)] 

(1.24) 

where F2{x,Q^) and R{x,Q^) (the ratio of longitudinal to transverse virtual 
photon cross-sections) are taken from parametrizations of unpolarized scat- 
tering data. Note that these induce a dependence in gi even if Ai is 
independent of . 

The possible reliability of the assumption that Ai is independent of can 
be explored using pertuiihat±Ye QCD models for gi{x,Q^), and several such 
studies have been made cJ. Leading-order analyses indicated a small 
dependence in Ai{x, Q^) (see eg. Ref. [§), but higher-order analyses question 
this (see eg. Refs. [|9|, ^],[^). The amount of any such dependence 
is sensitive to the polarization of the gluon density AG{x,Q^), and it may 
soon be possible, to use data to constrain this, though we do not believe this 
is yet reliable E3 (see, however, Ref. jsj]). It does not appear that such a 

dependence in g\'^{x,Q'^) would have impact on the integrals T\''^{Q'^) 
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outside the noted statistical and systematic errors. However, it could become 
an important effect in the future, and both theorists and experimentalists 
should keep their eyes open. 

The old-fashioned assumption of Regge behaviour at \ow-x also needs to 
be checked carefully. The leading-order perturbative QCD evolution equations 
for thejion-singlet part of the helicity distributions, Aq^si^jQ"^), lead us to 
expect L30 singular behaviour as a; ^ 0, so that 

A(jjvs(a;, Q^) ~ Cats exp(Ajvso- + Bns Inp - 7: In cr) (1-25) 

P 2 

where Ans: Bns and Cats are some constants and 



■»?'»^ ^^m- 

and we might expect by analogy with the unpolarized structure functions that 
xo ^ 0.1, Qg ~ 1 GeV and the leading-order QCD scale parameter A ~ 0.25 
GeV, with 

4\/2 4 

where Nf 3 in the range of current experimental interest (see also 
Ref. iQ). In principle Eq. (1.25) can be applied directly to the low- a; be- 
havior of the integrand of the Bjorken sum rule: gi{x,Q^) — gi{x,Q^) = 
i[Au(x,(3^) ~ Ad(a;, Q^)], as well as to the other nonsinglet combination 
Au{x, Q^) + Ad{x, Q2) _ 2As(a;, Q^) that also contributes to 5i'"(a;, Q^). The 
flavour-singlet combination of structure functions has a more complicated low- 
X behaviour, which could be important for the extraction of the Aq. It is not 
clear whether the behaviour in equation (1.25) is relevant to the data presently 
available: one SMC data point may be in its region of applicability, and could 
in principle be used to normalize the perturbative QCD formula, serving as a 
basis for extrapolating the integrals to x = 0. In practice, it does not seem 
at present that this would have a significant effect on the evaluation of the 
Bjorken sum rule. 

The analysis of the polarized structure function data has often assumed 
that the transverse polarization asymmetry 

^ da^^ — da^^ 
da^^ + da^^ 

is negligible. This is related to the spin-flip photon absorption asymmetry 

A, = (1.29) 

0'l/2 + Cr3/2 
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and the longitudinal Ai asymmetry (1.5) through the relation: 

A^= d (A2-7(1-|)Ai) . (1.30) 

were ai /2 and (T3 /2 are the virtual photon- nucleon absorption cross sections for 
total helicity 1/2 and 3/2, respectively, axL arises from the helicity spin-flip 
amplitude in forward photon-nucleon Compton scattering, 7 = 2AIx/ y/Q^ , 
and y — v / Eieptom where v is the energy transfer in the laboratory frame. 
The coefficient d is related to the virtual photon depolarization factor D by 



,.0^_ (1,31, 

The asymmetries Ai and A2 are subject to the following positivity conditions^ 

|^i|<l, 1^2! <Vi?. (1.32) 
and are related to the structure functions gi,2 by 

Ai = ^{gi - 7^52) , A2 = ^(51 + 92) , (1.33) 

where Fi = ^2(1 + -/'^)/2x{l + R) is the spin-independent^tmcture function. 
Recently, data on (72 have become available for the first timeocJ. They indicate 
that it is considerably smaller than the positivity bound in equation ( 1.32 ), 
and is very close to the leading twist formula of Ref . : 

gr{x,Q') = -g^{x,Q')+ J gi{t,Q')j , (1.34) 

e data on 172 are also compatible with the Burkhardt-Cottingham sum rule 

1 





g2{x,Q^)dx = 0, (1.35) 



which has been verified to leading order in perturbative QCDi§. The exper- 
imental errors are still considerable, in particular because the low- a; behaviour 
of g2 is less well understood than that of gi. However, the data already tell us 
that the uncertainty in 52 is not significant for the evaluations of the F^'" . 
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1-4 Higher orders in QCD perturbation theory: 



The perturbation series in QCD is expected to be asymptotic with rapidly 
growing coefficients: 

S{x) = y^Cnx", x= — ,Cn^n\K'^n^ (1.36) 

n=0 



for some coefficients iC, 7l3'E2I. This type of behaviour is associated with the 
presence of the renormalon singularities, as we shall discuss shortly. Such series 
are often evaluated approximately by calculating up to the "optimal" order, 
implicitly defined by 

|cn„,.a."-*|<|c„„,,+ia:"-*+i| (1.37) 

and assuming an error of the same order of magnitude as c„^pj jf-opt ^ ^j^g 
question arises whether one can approach or even surpass this accuracy without 
calculating all the terms up to order riopt —This possibility has been studied 
using the_effective charge (ECH) approach lIIlj and using commensurate scale 



relationsESo. Tn this section, we discuss the use of Pade approximants (PA's) 
for this purposeEJ'EJ. ■— n— ■ 

Pade approximantsEllS are rational functions chosen to equal the pertur- 
bative series to the order calculated: 

_ Op + aiX + ... + fljyX^ 

V^l^'^- l + 6ia:+... + 5MX*^ ■ 

(1.38) 

[N/M] = 5' + 0(2;^+*^+i) 



Under certain circumstances, an expansion of the PA in equation (1.3S) pro- 
vides a good estimate, c^^jv/+i' the Pade Approximant Prediction (PAP), for 
the next coefficient cp^^M+i in the perturbative seriesH. For example, we have 
demonstrated that if 

6„^^^-l^l, (1.39) 

as is the case for any series dominated by a finite number of renormalon sin- 
gularities, then (5[jv/jif] defined by 

A - - CN+M+l . , 

0[N/M] = (1-4U) 

CN+M+1 
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has the foUowing asymptotic behaviour 



^[N/M] 



where L = N + M + aM 



(1.41) 



and where a is a number of order 1 that depends on the series under consider- 
ation. This prediction agrees very weU with the known errors in the PAP's for 
the QCjD. vacuum polarization D function calculated in the large Nf approxi- 
mationE^I, as seen in Fig. 1.2a. 
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Figure 1.2: Relative errors in the [N/M] Fade approximants (a) to the |QpD vacuum 
polarization D-function, evaluated to all orders in the large- A'^j approximationEj (the rate of 
convergence agrees with expectations for a series with a discrete set of Borel poles), and (b) 
to the Borel transform of the D-function series, where the cn nverg -ence i s part icularly striking. 
The straight lines correspond to the error formulae, eqs. (1.41) and (|l.5l[), respectively. 



Large- A'^y^ calculations of the perturbative corrections to the Bjorken sum 
ruleEiJ indicate the presence of only a finite number of renormalon singularities, 
so that the PAP's should be accurate. Using the known terms in equation 



(1.16), the [1/2] and [2/1] PAP's yield the following estimates for the fourth- 
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order coefficient 



cf[^p^]«-lll ([1/2] PA) 



(1.42) 



^4 [PA] 



-114 ([2/1] PA) 



and the error estimator in equation (f .4f ) witli a = 1 yields 



'^[1/2] — -1/8; ^[2/1] 
These results can be combined to obtain 



4/4 



Bj 
"^4 [PA] 



-111 
1 + '^[1/2] 



-114 



1 



'[2/1] 



-139 



which is very close to the ECH estimate E 



-130 



(1.43) 



(1.44) 



(1.45) 



A second application of PA's is to "sum" the full perturbative series. The 
latter is ambiguous if the perturbative series possesses an infrared renormalon 
singularity, i.e. a divergence of the form in equation (1.36) with > 0. 
Consider the following toy example: 







1 - xt 



X 



-y/x 



1 



y 



dy 



(1.46) 



which exhibits an infrared renormalon pole at y = 1. One possible way to 
define the ambiguous integral on the right hand side of equation ( 1.46 ) is via 
the Cauchy principle value prescriptiontd. We see in Fig. 1.3 that the errors in 



the Pade "Sums" (PS's) [N/M]{x) are smaller than the truncated perturbative 
series 



N+M 
n=0 



(1.47) 



when n < Uopt, which is 5 in this example. You will notice in Fig. 1.3 that 



the errors in the PS's become unstable for large n: this is because of nearby 
poles in the denominator of equation (1.38) which are not important for small 
n. Also shown in Fig. 1.3 as "combined method" is a systematic approach to 
treating tliese poles and optimizing the PS's for large n, which is described 
elsewhere H 
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Figure 1.3: The relative errors between partial sums of the series S{x) = T,n\x" and the 
Cauchy principal value of the series (solid line) is compared with the relative errors of Pade 
Sums (dotted line). We see that the relative errors of the Pade Sums are smaller than 
those of the partial sums in low orders, fluctuate in an intermediate regime, and are again 
more accurate than the partial sums in higher orders. The fluctuations are associated with 
nearby poles in the Pade Sums, that may be treated by the "combined method" mentioned 
in the text, shown as the dashed line. 



Evidence that the PS's for the Bjorken sum rule provide a good estimate 
of the perturbative correction factor in equation (1.16) is provided by the 
study of the renormahzation scale dependence. We see in Fig. 1.4 that the 
renormalization scale dependence of the [2/2] PS is much smaller than that of 
the [2/1] and [1/2] PS's, which is in turn much smallp; than that of the naive 
perturbation series evaluated to third order. We recallo tha t the full correction 
factor should be scale- independent, and interpret Fig. lA as indicating that 
the PS's may be very close to the true result. 

Fig. 1.5 shows the estimates of the perturbative QCD correction to the 



Bjorken sum rule obtained in various different approximations, including third- 
order perturbation theory, fourth-order perturbation theory estimated using 
the ECH technique and the [2/1], [1/2] and [2/2] PS's. We interpret the latter 
as the best estimator, and take the difference between it and the [2/1] and [1/2] 
PS's as a theoretical uncertainty. Also shown in Fig. 1.5 is the experimental 
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Figure 1.4: The scale depende nce o f ai,(3GeV^) obtained from a fixed value f{x) = 
Co/qa) X 0.164 = 0.783, (c/. eq. ( |l.54| ) ) , for Q/2 < ^ < 2Q, using the naive third- and 
fourth-order perturbative series and the [1/2] and [2/2] PS's. 



error on this quantity, as extracted from the combined analysis of the available 
experimental data discussed in the next section. 

More information can be extracted by considering PA's in the Borel plane. 
The Borel transform of a perturbative series is defined by 

oc 

S'(x) = ^ CnX" 



Borel 
n=0 



~S{y)^Y.~^^y- : 5„ = £^ (1) (1.48) 



n=0 



n+l 



where /3o — (33 — 2Nf)/3. A discrete set of renormalon singularities would 
show up as a set of finite-order poles in this plane 



14 



1.0 



3rd order P.T. 



f(x) 



0.9 



O.f 



4th order P.T. - ECH 



[1/2] & [2/1] Pad 

[2/2] Pade - ECH 

0.7 ''''''''''''''' 
0.000 0.025 0.050 0.075 0.100 0.125 

x=a„/ 7T 




Figure 1.5: Different approximations to the Bjorken sum rule correction factor f{x), 
third-order and fourth-order perturbation theory, [1/2], [2/1] and [2/2] Pade Sums are 
compared. Also shown as a vertical error bar is the val ue oi f{x) we extract from the 
available polarized structure data (1.54). 



The PA's in equation (1.38) are clearly well suited to find the locations yu and 
the residues Vk of such poles. In the case of a perturbative series dominated 
by a finite set of L renormalon singularities, a sufficiently high-order PA will 
be exact 

[M/N]{y) = S{y) : iov M + N > Lq (1.50) 

for some Lq oc L. Generically, in any case where the quantity analogous to 
eq. ( 1.3£ ), e„ ~ the error analogous to ( 1.4C]| ) is given by 

(M!)2 



J[M/M] 



'2M 



(1.51) 



This prediction of very rapid convergence is confirmed in Fig. 1.2b in the 
case of the QCD vacuum polarization D function evaluated in the large Nf 
limit Lj, which has an infinite number of renormalon poles. 



The power of the application of PA's in the Borel plane is shown in Fig. 1.6 



where we see that the [2/1] PA of the Borel-transformed Bjorken series has 
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Figure 1.6: The locations and residues of poles in the [2/1] PA and in rational-function 
fits to the Borel transform of the first four terms in the perturbation series for the 
Bjorken sum rule. We note that the location of the lowest-lying infrared renormalon pole 
is estimated accurately by Pade Approximants in the Borel plane, and that its residue is 

stable in the different fits. 



a pole at y — 1.05 . The agreement with the expectation of a first infrared 



renormalon pole at y = 1 is striking. Fig. 1.6 also shows other rational-function 
approximations to the Borel-transformed Bjorken series. We see evidence that 
this is dominated by a strong infrared renormalon pole at y = 1, that there 
may be a weak ultraviolet renormalon pole at y = — 1 and possibly another 
pole at y = 2. 

The ambiguity in the definition of the perturbative Bjorken series asso- 
ciated with the y — 1 renprmalon singularity corresponds to a possible l/Q^ 
correction of magnitudeE3 



A(r5'-r^) =±i^ 0.98 7r^ (1.52) 



It is expected that the QCD correction to the Bjorken sum rule should include 
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a higher-twist correction of similar form with magnitude 



0.02 ±0.01 



(1.53) 



The perturbative ambiguity in equation ( 1.52| ) is cancelled by a corresponding 
ambiguity in the definition of the higher-twist contribution. In the next section 
we will treat equation (1.53) as a correction (with error) to be applied to the 
perturbative QCD factor shown in Fig. 
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We have also compared PA's to the predictions of commensurate scale rela- 
tions within the framework of Ref. |Q . The predictions of the two approaches 
are numerically very similar, and we give formal reasons in Ref. | |63[ | why we 
believe that this should be so. However, we shall not use commensurate scale 
relations in the data analysis of the next section. 



1.5 Numerical analysis of the Bjorken sum rule 

Table I shows the data on ±h£_inlfigrals F^'" currently available from experi- 
ments at CERN and SLACEaH . We do not attempt to correct these num- 
bers for any of the effects discussed in section 1.4, such as the Q^-dependence 
of the asymmetry ^i, the extrapolation to low- a;, or the transverse polarization 
asymmetry. We do not believe that any of these effects will change any of the 
data outside their quoted errors. We choose to evaluate the Bjorken sum rule 
at = 3 GeV^, which requires rescaling all the data as described in Ref. p6| . 

Table I 

pp.",d currently available from experiments at CERN and SLAC. 



experiment 


target 


Ti 


E142 


n 


-0.045 ±0.009 


E143 


P 


0.124 ± 0.011 


E143 


d 


0.041 ± 0.005 


SMC 


d 


0.023 ±0.025 


SMC ('94) 


d 


0.030 ±0.011 


SMC 


P 


0.122 ±0.016 


EMC 


P 


0.112 ±0.018 



All experimental data have been evolved to = 3 GeV^. 



17 



The following is the combined result that wc find for the Bjorken sum rule: 



r?(3GeV2) - V^iSGeV^) = 0.164 ± 0.011 



(1.54) 



which is indicated by a vertical error bar in the lower left corner of Fig. 1.5 
Comparing this value with the [2/2] PS estimate also shown there, we find 



(1.55) 



which becomes 



a,(M|)=0.119l°:°S^ ± 



(1.56) 



when we run as up to using the th ree- lo op re normalization group equation. 
The errors quoted in equations (|l.55| ) and (L5£) are purely experimental, and 
the second ± sign in equation ( |l.5(j ) indicates that further theoretical system- 
atic errors must be estimated. Those we have evaluated include that associated 



with the renormalization scale dependence shown in Fig. 1.4 (±0.002), the dif- 
ference between the [2/2] and [2/1], [1/2] PS's (±0.002), and the correction due 
to the higher-twist estimate in equation ( 1.53 ) (—0.003 ± 0.002), whereas the 
error in the running of as is fomid to be negligible. Combining these estimates 
with equation ( 1.56| ), we findE3 



a,(M|) = 0.116: 



-0.003 
-0.005 



±0.003 



(1.57) 



The stability of this result is indicated in Fig. 1.7, where we exhibit the 
values of as{M^) obtained using different orders of perturbation theory, com- 
pared with our result ( 1.57 ) obtained using the [2/2] PS. Also indicated is the 
shift induced by the higher-twist correction, which lies within our error bars. 

As can by seen from the compilation in Fig. |l.8| , our central value for 
as(Mf ) is quite compatiKb with other determinations and the world average, 
which is quoted to be BOO 



as{M^) = 0.117 ±0.005 



(1.58) 



Indeed, the error quoted in equation (1.57) is quite competitive with the most 
precise determinations of as{M'^) that are available. Moreover, we note that 
plenty of precise new data will soon be available from the SMC experiment at 
CERKL the E154 and E155 experiments at SLAC, and the HERMES experi- 
mentO at DESY. In the longer run, experiments with a polarized proton beam 
at HERA will provide valuable information on the behaviour of gi at low-x, 
as well as on its Q^-dependence at fixed x. 
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Figure 1.7: Values of (M^)_pbtained using different orders of perturbation theory, 
compared with our result (L.56), obta ined using the [2/2] PS. The size of the shift 
induced by higher-twist correction is ( |l.53| ) indicated by a downward arrow to the right of 

theT2/2] point. 
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from different sources (adapted from 



19 



1.6 Decomposition of the Nucleon Spin 

So far, we have only discussed the combination T\ — F" which enters in the 
Bjorken sum rule. The individual F^'" can be used as in equation (1.13), 
though not forgetting the perturbative QCD corrections in equation ( 1.17 ), to 



extract the individual Ag. As is seen in Fig. 1.9, the different experiments on 
both proton and neutron targets are all highly consistent, once the perturbative 
QCD corrections are taken into account. 



0.4 



0.2 



0.0 



-0.2 



average: 
0.27 ± 0.04 

= 2.0 



naive 
part on 
model 



0(a, ) 0{al ) 0{al ) 0{at ) 
X E142 + E143-P • E143-d O SMC-d(92) SMC-d(94) □ SMC-p W EMC 



Figure 1.9: The values of AE{Q^=3 GeV^) extracted from each experiment, plotted as 
functions of the increasing order of QCD perturbation theory used in obtaining AE from 
the data (from Ref. |66 , updated with most recent • 



Some time ago, it appeared as if the neutron data from E142 might be at 
variance with the other data points. However, this is no longer the case if all 
the higher-order corrections in equation (1.17) are taken into account, and the 
latest evaluations Il3 of the E142 data indicate a different preliminary value of 
F", as seen in Table F Making a global fit, we find 
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Au = 0.82 ± 0.03 ±... 

Ad = -0.44 ± 0.03 ±... (1.59) 
As = -0.11±0.03±... 



and 

AI]== Au + Ad + As = 0.27±0.04±... (1.60) 

where the second ± sign indicates that further theoretical and systematic er- 
rors remain to be assigned. These include higher-twist effects, errors in the 
extrapolation to low- a; which is more complicated than for the nonsinglet com- 
bination of structure functions appearing in the Bjorken integrand, the possible 
Q^-dependence of Ai, etc.. We believe that th ese er r ors m ay combine to be 
comparable with the errors quoted in equations (1.59),( i.60| ), but prefer not to 



quote definitive ranges for the Aq until all these errors are controlled as well 
as those appearing in the Bjorken sum rule. 

The|_SMC Collaboration has presented at this school a global analysis of 
the Aq llJ, in which they combine SLAC and CERN measurements of the 
structure functions to arrive at new estimates of the F]"'", which differ from 
the values we have used (see Table I) principally because of their treatment of 
the \ow-x data. An analysis of a subset of the the available data has also been 
presented in Ref. ^J]. Both these analyses differ from ours in the treatment 
of higher-order perturbative QCD effects, but agree within the stated errors. 

One may also get a feeling for the expected range of AE and As by plot- 
ting the results for these two observables extracted from each of the existing 



experiments, as shown in Fig. 1.10 



1.7 Outlook 

In this lecture we have concentrated on the phenomenological analysis of the 
data on polarized structure functions presently available. As we have seen, 
these tell a remarkably consistent story, once higher-order QCD corrections 
are included. We have not addressed in great detail here the theoretical in- 
terpretation of the data, nor their spin-offs in hadron physics and elsewhere, 
nor possible future developments in this field. In fact, these measurements 
provide valuable insights into important issues in ncitLjperturbative QCD, such 
as the role of chiraL sjonmetry in nucleon structureEj, the axial anomaly and 
the [/(p,ip£pbleml2l~E3, and the relationship between current and constituent 
quarks which are provoking lively theoretical debates (see Ref. for a 
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Figure 1.10: The values of AS and As extracted from each experiment, plotted 
against each other. All data have been evolved to common q2 = 3 GeV^. The clear 
linear correlation between AS and __As _tcaults from the linear relations 

(g),(i.ic),(i.i7). 



recent application to the pion structure). Some of these issues are reviewed in 
Lecture II. 

The polarized structure function data support previous indications from 
the TT-nucleon tr-term and elsewhere that strange quarks in the nucleon wave 
function cannot be neglected, with interesting implications for the analysis 
recent data from LEAR on production in proton-antiproton annihilation 
, as also discussed in Lecture II. Among other spin-offs, we recall that the 
axial-current matrix elements extracted from polarization data determine scat- 
tering matrix elements foL candidate dark_matter particles such as the lightest 
supersymmetric particleEj and the axiontj, as discussed in Lecture III. 

In the future, we look forward to the completion of the SMC programme 
and its possible HMC successor at CERN, the E154 and E155 experiments 
at SLAC, data from the HERMES experiment at HERA, the polarized pro- 
ton programme at RHIC, and possible polarized electrons and protons in the 
HERA ring. The tasks of these experiments will include the determination 
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of the gluon contribution to the nucleon spin, the elucidation of the Q,^ de- 
pendence of Ai , and the \ow-x behaviour of gi . These will continue to fuel 
activity in this interesting field for the foreseeable future, which will lead us to 
a deeper understanding of the nucleon, an object we thought we knew so well, 
but which reveals a new face when it spins. 
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2 Chiral Solitons and Strangeness in the Nucleon 

In this lecture, we review possible interpretations of the polarized structure 
fiinct inn data, with particular emphasis on the chiral soliton (skyrmion) pro- 
posal EEl, but with comments on other possibilities EJl3. We also discuss in- 
dications from LEAR for a pattern of violations of the Okubo-Zweig-Iizuka 
(OZI) rule in (p production in pp annihilation, motivating a phenomenological 
model for ss component in the nucleon wave function, whose possible tests we 
also review. 



2. 1 Chiral soliton interpretation of EMC spin effect 

In this section, we will describe a possible interpretation of EMC spin effect in 
terms of a chiral soliton model of hadrons. This ha^i the double interest of being, 
as far as we know, the only model which explainsc3 the small experimental value 
of AS and can also be derived from the underlying QCD theory, as we now 
describe (see Ref. Q for a general introduction to QCD and its symmetries). 
The classical QCD Lagrangian 

C-QCD = Cym ~ ^{qLpqL + QRpqR) - ^mgiqLqR + qRqh) (2.1) 

q=l q=l 

is invariant under the following global chiral transformations 

qL LqL and qa RqR; L, R e U{Nf) (2.2) 

which constitute the symmetry group SU{Nf)L x SU{Nf)fi^ x U{l)v x U{1)a, 
where V = L + R and A = R—L. As will be discussed later, the U{\)a part of 
the chiral symmetry is broken by quantum effects, and we concentrate for now 
on the rest of the chiral symmetry group. It is believed that non-perturbative 
QCD dynamics break this chiral symmetry spontaneously to the subgroup 
SU{Nf)v- In agreement with the Goldstone theorem, this spontaneous chiral 
symmetry breaking is accompanied by the appearance of N'j — 1 Goldstone 
bosons, which are the tt^'", K^''^, K^, rjs in the case of SU{3) with the u, 
d, and s quarks. In the limit of small u,d,s quark masses, where SU{2i)v 
symmetry becomes exact, the interactions of these light pseudo-scalar bosons 
are described by the following effective Lagrangian: 

|Tr(a,f/a-C/t) + J-, 



Cps - (^Tr(9^[/a^C/t) + — Tr[(a^C/)[/t,a''C/C/t]' 



(2.3) 
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where J7 is a matrix which represents the pseudoscalar Goldstone boson fields 



a= 1. 



U = exp 



(2.4) 



The first term in ( |2.3| ) provides the P-wave rnr scattering lengths, as well 
as the canonical kinetic term for the boson fields. The second term contains 
higher order in the field derivatives, such as may be generated by vector-meson 
exchange, or by quark loops. The third term is related to the light quark 
masses, and yields 



A, 



QCD ruq . 



(2.5) 



When Nf > 3 there is an additional Wess-Zumino term Cwz which does not 
concern us here. 

The above effective Lagrangian (2^) is believed to clescribe accurately 
QCD dynamics at energies E <^ ^qcd for the light quarksLj: 



rUuM ^ 10 MeV, m,, < 150 MeV 



(2.6) 



It may be regarded as the first step in a systematic bosonization of QCD, i.e., 
the expression of QCD dynamics in terms of meson fields. This bosonization 
can be derived formally in the l/N^ expansion, which is often treated in the 
lowest-order approximation, though this is not necessary. In the l/Nf. expan- 
sion, one considers a replacement SU{S)c — > SU{Nc), in which 



127r 



(2.7) 



{nNc-2Nf) lnQ2/A2 
We see explicitly that for large Nc the gauge coupling strength ^ l/A^^. In 






Figure 2.1: Planar diagrams for mesons in the large- A^c limit. The solid lines in the right- 
hand part represent colour lines, and the dashed line flavour. 

this expansion, QCD is regarded as a theory of an infinite number of mesons, 
whose dynamics is described in leading order by planar diagrams such as those 
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in Fig. 2.1, It is easy to verify that these are aU of the same order in the hmit 
Nc oo, g^Nc fixed. Higher orders in the 1/Nc expansion may be treated 
systematicaUy using diagrams of higher topological order. In two dimensions 
this programme hasJpeen realized explicitly; indeed QCD2 has been "solved" 
for any value of Nc c3. 

All well and good, but where are the bajryons? The answer is that they are 



solitons in this apparently bosonic theory Hti. The ener gy o f any field con- 
figuration described by the effective Lagrangian equation ( |2.3| ) can be written 
as 

E = J d^xi^^TT0U-dU^) + ^TT[{d,U)Ul{d,U)U^^ 

This expression is finite if the following condition is satisfied: 

U{x) -^1 as |f I cx) (2.9) 

Elementary topology tells us that such field configurations are classified by the 
group U3{SU{Nf)) = Z, which characterizes all the possible ways of mapping 
the 3-sphere into the chiral symmetry group. These different field configura- 
tions may be labeled by the quantity 

B^^^vk J d^'xTvid.UU^djU^dkUU^) (2.10) 

which a topologist would identify as the winding number of the configuration. 
The quantity B is related in the normal way to a conserved current: 

B = J d^xJo : 9^J^ = (2.11) 

As discussed extensively elsewhere @@, B can be identified with baryon num- 
ber. In the case of a spherically-symmetric soliton configuration: 

U(x) ^ cos9(r) +i^—^siiie{r) : r = \x\ (2.12) 
r 

one has: 

B = - [9(0) - 0{oo)] (2.13) 

TT 

which is unity if 9{0) — vr and ^(oo) = 0, which i^te case for the lowest-energy 
soliton with B = 1. Also as discussed elsewhereEaa one can demonstrate that 
this B = 1 soliton is a fermion. 
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The above discussion may sound rather high-falutin', but the construction 
nfsuch chiral solitons has been carried out explicitly in QCD in two dimensions 
O as we will discuss shortly. 

An essential complication in the treatment of these chiral solitons is the 
fact that there are many equivalent soliton configurations: if Uq is a solution, 
then so is 

Va = VUoV~^ (2.14) 

where V is an arbitrary spatially-constant SU{Nf) matrix. It is essential to 
take account of this when one quantizes the theory. This is done by parametriz- 
ing V in terms of its internal variables, the collective coordinates, and consid- 
ering soliton wave functions x(^) in the space of such collective coordinates: 



\N) = J dVx{V)\V) (2.15) 
At the quantum level, any state |A^) must have definite isospin / and spin J, 




Figure 2.2: A one-dimensional analogue of the collective coordinate V: particle constrained 
to move on a circular ring. Classical ground state corresponds to a particle at rest at some 
fixed angle 9. In quantum mechanics this is no longer true and we must have an eigenstate 
of the angular momentum operator Lg = —i-^. 

and this requires a superposition in the internal V space. The simplest way 
to see this is to consider the analogous case of a point particle located on a 
ring: Uq = e'^ in this case. Classically, the particle may sit at any fixed value 
of the angle 9. However, at the quantum level the only consistent state is an 
eigenstate of angular momentum, i.e., the point particle must rotate around the 
ring with an angular velocity ui, which is described hy V = e*"* ! Exactly the 
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same phenomenon occurs with the SU^Nf) chiral sohton: the quantum wave 
function must correspond to a rotation in the internal V space. This rotation 
Lu is slow (~ l/-^c) in the large- A'c limit, in which a semi-classical treatment 
of the quantum state is sufficient, but one need not in principle restrict oneself 
to this assumption. In thp limit; of large Nc, the physical parameters of the 
soliton behave as follows oE3 : 

m^-N,, Rn-N^, I^N,, J^i^iVO (2.16) 

where i?jv is the nuclcon radius, and / its moment of inertia. The correspond- 
ing picture of the baryon that emerges is one of a rotating collective state of 
pseudoscalar mesons. The rotation velocity to — J / 1 is slow in the limit of 
large Nc^ but the underlying physical picture is valid in principle for any value 
on Nc- 

As already mentioned, the above picture has been realized explicitly in 
QCD in two space-time dimensions for arbitrary values of Nc and A^/O. The 



solitons have precisely the form of equation (2.14), with the matrix [/o taking 
the form 



/I 



C/n = 



exp 



(2.17) 



where 



. / N /4A^c _i 
00 (x) = -\/— — tan 



exp \ / — mx 



(2.18) 



The SU{Nf) rotation matrix V can be written in the form 



V = 



( ^1 \ 

V 



(2.19) 



V ZNf I 
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where the form of V does not concern us here, while the parameters 
Zi : i — 1, . . . , Nf play the role of the collective coordinates introduced ear- 
lier. They characterize the soliton solutions as follows: 

i^(z,,z*) ^ X Z2"^ X ... X zly X zr' X . . . z^^""/ : ^ (n, ~ m,) - 

i=l 

(2.20) 

The lowest-lying states lie in a 10 representation of SU (3) analogous to that 
containing the A and 57 baryons in four dimensions. Thus we indeed have an 
explicit realization of the soliton ideas discussed previously. In particular, it is 
possible in this model Jx) calculate the ratios of the matrix elements {B\qq\B) 
in these baryon statesEj. In the general case of Nf flavors and Nc colors, one 
obtains 

{msea)B = TT^-TT' (2-21) 

where {qq)sea refers to the non-valence quarks in the baryon B. One can also 
compute flavor content of valence quarks. Consider a baryon B containing k 
quarks of flavor v. The i;-flavor content of such a baryon is 



1 



Nf + Nc 



(2.22) 



For Nf =3, Nc^3 one finds 



{A++\uu, dd, ss\A++) (X 4,1,1 

(2.23) 

{fl^\uu, dd, ss\fl^) cx 1,1,4 

We see explicitly that the ss matrix element in the non-strange baryon A^"*" is 
not negligible, and neither are the uu, dd matrix elements in the triply-strange 
baryon Q~ . 

Unfortunately, it is not possible to explore axial-current matrix elements 
directly in this two-dimensional model, so we return to four dimensions. 

In the picture of spontaneously-broken chiral symmetry, the matrix ele- 
ments of axial currents are given in general by the PCAC hypothesis, namely 
that they are dominated by couplings of the pseudoscalar Goldstone boson 
with the same quantum numbers. We know from tt^ decay that 

{0\A^\7r) =iUpf, (2.24) 
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and similarly for the , and the same is believed to be true for other members 
of the pseudoscalar-meson octet. Matrix elements of the SU{3) octet axial 
currents are related to the couplings of the octet pseudoscalars: 

{X\A'\Y) ^ Ug^.^^XY ■ (2.25) 

The ninth axial current is special, because it is not conserved at the quantum 
levelB: 

d^Al cx g^F^^.F'"' (2.26) 

For this reason, the ninth pseudoscalar rf is onlv_a pseudo-Goldstone boson, 
whose mass would be non-zero even if mu^d,s = OlZI : 

m„„ cx 1/iVe (2.27) 

It should also be noted that the ninth pseudoscalar decouples from the other 
eight in the ILaiit of large N^. since its couplings to them proceed via interme- 
diate gluonsEj : for example, there is a coupling 

'C.,.o = ^Tr(a^C/9'^C/t)^ (2.28) 

As already commented, the chiral soliton exists because of the topology of 
SU{Nf), not that of U{1). This means that the baryon is to be regarded as a 



"lump" of the n/K/r/g, mesons, not the singlet 770. Moreover, equation (2.28) 
indicates that the 770 decouples from the baryon at leading order in l/Nc- 

Chiral soliton calculations are generally made in the relatively easy limits 
Nc 00, mu.d,s 0. One should beware of TVc-dependent predictions, which 
depend on the truncation of the a priori infinite-dimensional meson theory. 
Examples of such bad iVc-dependent quantities include 

: > gA , /tt , g-RNN , g-KNA , ■ ■ • (2.29) 

On the other hand, calculations of A'c-independent quantities such as 

{r'^)N , Mp//^n, filgA, gTTNA/gTTNN (O 'in\ 

(19%) (2%) (3%) (1%) ^ """^ 

may nnt be so unreliable, as indicated by the accuracy of the predictions listed 
above l3. Another example of a quantity which may be reliably estimated in 



the SU{3)f limit isl; 



(p\-ss\p) 7 ^2.31) 



{p\uu + dd + ss\p) 30 
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and this number is indeed consistent with determinations of the 7r-nueleon a 

strated successes in calculations 



ornrtr. 



LOC 
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term. The chiral soliton model has also 
of meson-baryon scattering phase shifts 

The chiral soliton mOjdfil offers three ways of seeing that |p) ~ 0, i.e. 

that Aw + Ad + As ~ cEl. The brute- force method is by direct calculation 
using the following standard representation of Al^{x) in the soliton state 

J d^xAl{x) oc Tr [XaV-^X,V] (2.32) 

it is easy to see that the matrix clement of vanishes because Tr(Aa) = 0. 

The second way to see the vanishing of (p|ylj^|p) is to consider the Goldberger- 
Treiman relation for the baryonic matrix element of the ninth axial current. As 
discussed above, the PCAC hypothesis relates this to the baryonic coupling of 
the ninth pseudoscalar meson rjQ. To see that this coupling vanishes, it suffices 
to make another trivial trace over SU{Nf) indices: 

g^,NN (X Tr (A,J7(f)) = Tr (^X.VUoV-'^ (2.33) 

This vanishing result can be understood at a more basic level by remembering 
that the baryonic soliton consists of tt, K and rjs mesons alone in the large- iVc 
limit, and that the rjQ meson decouples from the others in this same limit. 

The third way of seeing that (p|A°|p) ~ is to consider the underlying 
physics of the soliton spin, by analogy with the particle on a ring discussed 
earlier. At the classical level, the baryon mass is given simply by 

M = y" d^x&oo{x) (2.34) 

where the energy- momentum tensor is derived from a static solution Uo{x). 
However, we recall that this is not an eigenstate of spin or isospin: these rotate 
with an angular velocity lu ^ l/Nc- The resulting expression for the soliton 
angular momentum is 

Ji = / (fxeijkXj 9ofc(x) = / (2.35) 



where 

T = 

3 

is its moment of inertia. Including the corresponding rotational energy, equa- 



1=1 I d^xOoor^ (2.36) 



tion ( 2.34 ) is modified to become 
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where the quantization of angular momentum imposes 



J(J + 1) = ^ for N, ^ for A (2.38) 

We see exphcitly from this construction that all the nucleon angular momentum 
is orbital in origin. If we consider the angular momentum sum rule 

^ = ^E^9 + AG + i. (2.39) 

this argument tells us that 

AI] = ^Ag = 0, = i (2.40) 

q 

Simple niDdels which extend the chiral soliton model to include gluons suggest 
also thatB 

AG = (2.41) 

These results have been derived in the limit of large Nc and small niu^d.s- At- 
tempts have been made to calculate corrections to this double limit, but these 
are incomplete so far. Nevertheless, nai've guess and incomplete calculations 



suggest that the result equation ( 2.40 ) might be accurate to within 30% or so. 
Since this range includes the present experimental value discussed in Lecture 
I, we consider this an encouragixLg success for the chiral soliton model. 

An alternative approachcil c3 to interpreting the polarized structure func- 
tion data is based on the axial U{\) anomaly, which should be taken into 
account in computing matrix elements of the singlet axial current. The defini- 
tion of the polarized quark distribution lS.q{x, Q^) is ambiguous at the 1-loop 
level. One possible interpretation replaces 

Aq Aq= Aq- -^AG (2.42) 
27r 

where the second term on the right-hand side is derived from the gluon con- 
tribution AG to the proton spin. In this interpretation all the previous de- 
terminations of Aq should be rephrased as determinations of the Aq, and it 
is possible, in principle, to resurrect the original quark model assumption by 
postulating that 

As = As - ^AG 7^ : As = 0, AG ^ (2.43) 
2n 
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In order for this mechanism to work, the value of AG must be quite large 



AG - 2 



(2.44) 



at ~ 3 GeV^ . This requires what might at first sight appear to be a rather 
bizarre decomposition of the proton angular momentum: 



+ (.2) + (.-2) 



(2.45) 



However, it should be noted that a leading-order compensation between large 
values of AG and is natural in perturbative QCD. Nevertheless, this anoma- 
lous gluon mechanism does not really explain the EMC result, though it may 
accommodate it. It is not yet possible to test definitively this model. For 
one thing, the x-dependence of the possible AG contribution requires further 
specification. So far we are aware of just one experimental result which bears 
on AG: a Fermilab fixed-target experiment has searched foii^ asymmetry in 
multiple 7 production in polarized proton-proton collisions!] 
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Figure 2.3: The particle production asymmetry Ai^i^ measured in Ref. [L03|, compared 
with various phenomenological models labeled by the curves (I) to (VI), as functions of the 
effective mass M of the produced system. 
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As seen in Fig. 2.3, they see no significanL^ymmetry, which is in conflict 
with some polarized-gluon models, but not allE23. Unravelhng the polarization 
of the gluons in the proton is one of the primary objectives of the next round 
of polarization experiments such as the HMC and RHIC. 
J—. pLn addition to the soliton and anomaly interpretation, it was also suggested 
E3~Ea that a significant suppression of the QCD topological susceptibility might 
play a key role, which would modify the naive quark model predictions (see 
also the discussion in Lecture I). 

2.2 The Okubo-Zweig-Iizuka Rule 

The polarized structure function experiments are not alone in indicating that 
there may be strange quarks inside the proton wave function. Other indications 
come from experimental violations of the Okubo-Zweig-Iizuka (OZI) rule 1123, 
as discussed in the rest of this lecture. 

According to the OZI rule, the only strong-interaction processes allowed 
are those which can proceed via connected quark line diagrammes such as those 



shown in Fig. 2.4 




Figure 2.4: Connected quark line diagrams for p decay and the nnA coupling, using the naive 
quark model content of the particles, and hence corresponding to processes conventionally 

allowed by the OZI rule. 

This assumption must be supplemented by Ansatze for hadron wave func- 
tions, such as the naive quark model contents 

= \p)^\uud) (2.46) 

Indications from meson mass formulae and the observation of — > 37r decay 
are that the wave function is actually a mixture: 

10) =cos(5|ss) + sin 6 \uu + dd)/V2 (2.47) 

where (5 <C 1. Disconnected diagrams may be mediated by gluon exchange, as 
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b) 



c) 



Figure 2.5: Disconnected diagrams mediated by gluon exchange, relevant to (a) (j> and J/ip 
decays, (b) a departure from ideal mixing in the ip wave function, and (c) transitions between 

different charmonium states. 



shown in Fig. |2.5| , which are subject to dynamical suppression which depends 
on the process considered: 



94,31: 



10' 



10" 



^had 



10- 



10" 



(2.48) 



The OZI rule finds some justification in the large- A^c limit of QCD considered 
earlier. Simple Nc power-counting for the meson diagrams in Fig. ^.(\ shows 
that they are suppressed by 1/Nc and l/N^, respectively. 




Figure 2.6: Meson diagrams that are suppressed in the large- A^c limit. As in Fig. 2.1, the 
solid lines represent colour, and the dashed lines flavour. 



ilawever, the applicability of the OZI rule to baryons is more question- 



ableliilj. It is more difhcult to justify in the large Nc limit, because the baryon 



wave function contains 0{Nc) quarks. Moreover, it does not seem to work very 
well experimentally. Historically, one of the . fir ^ t ip dications of a discrepancy 
with the OZI rule was the 7r-nucleon cr-termt 
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= ^(™« + md){p\uu + dd\p) 



(2.49) 



Using the Gell-Mann-Okubo mass formula and the OZI assumption that {p\ss\p) 
one estimates 

S"'^ ~ 25 MeV (2.50) 
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to be compared with the experimental value of about 45 MeV, corresponding 
to 

{p\ss\p) 



0.2 



(2.51) 



{p\uu + dd + ss\p) 

The data on polarized structure functions discussed in Lecture I provide an- 
other example of an ss matrix element in the proton which is non-negligible. 

In addition to these static matrix elements, there have long been indica- 
tions of OZI "violation" in the production of the 4> and other supposedly pure 
ss mesons. For example, old bubble chamber data indicated that 



(7{pp — > (flTT^TT ) 

(y{pp — > (/)7r+7r~) 



(19 ±5) X 10" 



(2.52) 



which is much larger than what could be expected from an admixture of uu, 
dd components in the (f) wave function, which is estimated to be around (1 
4) X 10^'^ on the basis of mass formulae. The bubble chamber result eq. ( ^.52 ) 
corresponds to 



0.05 < 



V2A{_ 



pp 



X) 



A{pp ^uu + X)+ A{pp dd + X) 



< 0.22 



(2.53) 



In addition to this observation, early data indicated possible OZI violations in 
biT~ , pp — > /2(1520)7r+7r~, pp ppcj) . . . and Q.* Vbui decay. 



123 due to an ss 



pn 

Among the proposed interpretations was OZI "evasion" 
component in the proton wave function, which provides the possibility of a 
new class of connected quark line diagrams as seen in Fig. 2/7. Additional 
evidence for this hypothesis may come from the presence of a backward peak 
in the reaction pp — > K" which could be due to the connected quark line 
diagram shown in Fig. 2.8, which involves the intrinsic strange component of 
the proton wave function. 

An alternative explanation for "excess" (f> production that has been pro- 
posed is rescattering through intermediate K, K* states, but it has also been 
argued that any such effect would be too small, and subject to systematic 
cancellations. Yet another interpretation was proposed in terms of an exotic 
1 resonance C(1480) previously reported in the ^tt*^ mass spectrum in the 
reaction 7r~p Cn. However, it seems difficult to reconcile this interpretation 
with the final-state channel dependence seen more recently, and LEAR data 
do not confirm the existence of such a resonance. 



36 





Figure 2.8: New class of connected quark line diagram which could account for the appear- 
ance of a backward peak in the reaction pp K~ . 



It is also worth noticing that dispersion relation analyses of the proton 
vector ispsealac, form factor also found suggestions of a surprisingly large (j)pp 



couphngll23~ti3 : 



2 



0.211 to 0.276 (2.54) 
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2.3 New Data from LEAR 



The LEAR ring at CERN has provided large numbers of nuclcon-antinucleon 
annihilations at or close to rest, which have enabled the OZI rule to be tested in 
(f) production in association with many different final states X . The ASTERIX 
collaboration has measured the ratios of (t)X to loX production for the states 
X = 77,7], Lu, p and tttt, as reported in Table II. 



The ratios R - 



TABLE II 

4>X/ujX for production of the (p and uj - 



mesons in antinucleon annihilation at 




rest. The parameter Z of the OZI-rule violation is calculated for 5 = — = 3.7 , assuming 
identical phases of the 4> and uj production amplitudes. The data are given for annihilation 
in liquid hydrogen target (perceiptage of annihilation from P-wave is ~ 10 ^ 20%), gas target 
(~61% P-wave) and LX-triggerEJ (~86-91% P-wave). 



Initial states B.R.-IO* 



Final state 



B. ■ 10^ 



\Z\ (%) Comnye^s 




07 



0.17 ±0.04 250 ±i 



42 ±8 



liquid, 



5.5 ± 0.7 
1.9 ±0.5 
0.3 ± 0.3 



96± 15 



24 ±2 



9.0 ± 1.1 
14.8 ± 1.1 



83 ±25 
133 ± 26 
113 ±30 
110 ±15 



22 ±4 
29 ±3 
27 ±4 
26 ±2 



<t>ri 
<t>V 
4>r) 



0.9 ± 0.3 
0.37 ±0.09 
0.41 ±0.16 



6.0 ±2.0 



1.3±1.2 



3.4 ±0.8 
4.4 ± 1.2 



6.3 ± 1.6 
7.5 ±2.4 



1.4 ±1.0 
2.1 ± 1.2 



6.3 ± 2.3 
3.0 ± 1.1 
4.2 ± 1.4 



19 ±7 



7±4 



PX,X : 

<px,x ■■ 



■ Pj 1.2 ±0.6 

4.6 ± 0.9 
5.4 ± 1.0 

7.7 ± 1.7 



6.0 ±3.0 
7.0 ± 1.4 
7.9 ± 1.7 
11.0 ±3.0 



1.3 ±2.0 
1.9 ±0.8 

2.4 ± 1.0 
4.0 ±1.4 



Tp < 200 MeV/c 
* p > 400 MeV/c 



LX-trigger 



We see that there are large enhancements in S'-wave annihilations as compared 
to naive OZI rule, especially for the case X = tt, whereas there are little or no 
enhancements in P-wave annihilations for A" = tt, 77, p, tttt. The Crystal Barrel 
collaboration has measured the corresponding ratios for X = tt*^, ry, tt'^tt*', 7. 
As also seen in Table II, they confirm an enhancement over the OZI rule in 
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the case X — tt'^, see no big effect in the cases X ~ 77, tt^tt^, and find a very 
large enhancement in the case X = 7, which is about a hundred times the OZI 
prediction! The OBELfX collaboration also finds large enhancements in the 
cases X = TT*, as also seen in Table II. 




Figure 2.9: Recent experimental evidence from LEAR for breakdowns of the naiVe OZI rule, 
(a) A </) peak is seen clearly in annihilations at NTP, but is less prominent at a pressure of 5 
mbar, indicating (b) that 4>-k production is suppressed in P-wave annihilations. The amount 
of OZI "violation" depends on (c) the invariant mass of the system produced in association 
with the meson, and (d) the invariant momentum transfer to the 0. 

Some interesting features of the available data are shown in Fig. |2.9|. We see 



in Fig. 2.9a that </)7r production is much smaller in P-wave annihilations than in 
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5- wave annihilations, and Fig. 2.9b shows that the tttt spectrum also depends 



on the partial wave. Fig. |2.9| c shows the dependence of (f>X/ujX production 
ratios on the invariant mass of the system X , where we see progressively larger 
enhancements at smaller masses. Fig. 2.9d shows the corresponding ratios 
plotted versus the invariant momentum transfer where we see a large effect 
at large t. 

These very interesting data exhibit the following features which need to 
be explained or accommodated in any model of OZI violation: 

- Non-universal enhancement factors, which are strong for X ~ tt, 7, 
smaller for X — p,LO^ -kt:, and not apparent for X ~ rj. 

- Larger enhancements in proton-antiproton annihilation at or near rest 
than in higher-energy annihilations, or in tttt and pp scattering. 

- When the initial pp state is known, large enhancements are seen in S- 
wave annihilations into X ~ n, tttt, but not in P-wave annihilations. 

The next section introduces a modell£3 which accommodates and explains 
these key features, and also makes some predictions for possible future exper- 
iments. 



2.4 Model for a Polarized Strange Component in the Proton Wave Function 

We consider ii the likelihood that the Fock-space decomposition of the proton 
wave function contains an ss component which we parametrize as follows: 

|p >= x \uudX) + z \ uudssX) (2.55) 

X X 

where X denotes any combination of gluons and light qq pairs. As seen in 



Fig. 2.10, two important new classes of connected quark line diagrams become 



possible, the shake-out diagrams illustrated in Fig. 2.1Ca and the rearrange- 



ment diagrams illustrated in Fig. 2. Kb. A typical shake-out amplitude may 

be estimated as 

Aipp ^ ss -I- a:) ~ 2Re{xz*) P{ss), (2.56) 

where P(ss) is a projection factor which may depend on the final state, and 
perhaps also the initial state. A typical rearrangement amplitude may be 
estimated as 

Aipp ^ ss + X) \z\^ T{ss) (2.57) 

where T{ss) is another projection factor which is very likely to depend on the 
initial state: for example, naive quark models might suggest that (f> production 
would be enhanced in rearrangement of an initial 5- wave pp state. 
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Figure 2.10: "Shakcout" and "rearrangement" diagrams responsible for {pp — > ss + X) in 
presence of an ss component in the proton wave function. 



If we define the generic amplitude ratio 

AiA + B^uu + X)+ A{A + B^dd + X) ^ ' 

the measured production ratios are given by 



_ (7{(I)X) ( Z + tan(5 



Rx — Y ~ ( r ' """" ) X (phase-space ratio) (2.59) 
a{uoX) \\-Zi&n5J ^ > ^ > 

where b is the angle representing the departure from ideal mixing in the vector 
meson nonet, which is a measure of the "expected" deviation from the OZI 
rule. We see from equation ( 2.55| ) that one might expect generic shake-out 
diagrams to yield 

\Z\ = 2- ^2 (2.60) 

assuming similar projection factors P, and generic rearrangement diagrams to 
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yield 



(2.61) 



assuming similar projection factors T . Data on S'-wave tt production corre- 
spond to 

|Z(07r/w7r)| = 0.24 ± 0.02 (2.62) 

leading to the prediction 

0.01 < |zp < 0.19 (2.63) 

to the extent that the reactions are dominated by the shake-out and rearrange- 
ment mechanisms, respectively. These estimates are both compatible with data 
on K production in pp annihilation, which can be expected to contain a con- 
tribution from the higher Fock-space states in (2.55) as well as ss production 
in the final state: 



Yk - (4.74 ± 0.22)% ~ A[Re{xz*)f 



(2.64) 



U 
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Figure 2.11: Emission of a /C+ meson in a chiral model, which may explain the negative 
polarization of s quarks in the proton wave function. 

As we have seen in Lecture I, polarized structure function data indicate 
that the strange component of the proton wave function is polarized negatively: 



As < 0, and this feature should be incorporated in ( 2.55 ). Negative pol ariza- 
tion of s quarks is understandable in a chiral model, as seen in Fig. 2.11 
The emission of a if"*" flips the positive helicity of a m quark into negative 
polarization for an s quark. This argument does not provide immediately an 
indication on the possible polarization of an ss antiquark, but the simplest 
allowed possibility is 

|W)i 1 |ss)i,-i (2.65) 

where there is a relative orbital angular momentum |1,1) between the two 
indicated components. The internal state of the ss pair is motivated by the 
fact that condensation in the vacuum occurs in a ^Pq state. 
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Within this picture, one could expect that the rearrangement diagrams 

from a spin-triplet initial pp state would yield preferentially spin-triplet ss 
states such as the </>. One could also expect that 5*- wave pp annihilations would 
yield preferentially final states containing S-wave ss pairs. These arguments 
favour maximum (j) production in ^5*1 annihilations, and less enhancement in 
^5*0 and P-wave annihilations. A corollary would be the dilution of the <p 
enhancement at higher energies, where more partial waves contribute. Pre- 
dictions for final states containing rj mesons are more complicated, because 
more diagrams contribute since the rj wave function contains both qq and ss 
components. 

2.5 Future Tests of the Model 

The following are some possible tests of these ideas for OZI evasion that could 
be tested in future experiments. 

1. Production of the P-wavc quark- model state /'(1 525) could be enhanced 
in P-wave pp annihilations, at least to the extent that they receive con- 
tributions from rearrangement diagrams. This suggestion already seems 
to be compatible with preliminary LEAR data. 

2. A reduction in the enhancement of ^tt/wtt as the p momentum is in- 
creased, associated with the decreasing fraction of S-wavc annihilations 
- for example, the 5- wave fraction at Pp = 600 MeV is between 14 and 
20 %. This suggestion also seems to be compatible with data from the 
Crystal Barrel collaboration. 

3. There should be considerable spin-dependence in the cross section for 
pp (pcj), which should, for example, be higher from spin-triplet initial 
states. 

4. Both S- and P-waves may contribute to the (/>7r7r final state, and we 
would expect dominance by the initial state, which seems once more 
to be supported by recent data. 

5. We would expect strong spin correlations in pp — > K*K*, which should 
be dominated by the L = 0, S = 2 state. 

6. The large enhancement of cjrf/uj^ observed should be dominated by the 
^fo,i,2 initial states, rather than the state. 
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7. The angular distribution of the e~^e pairs produced by (j) decay in reac- 
tions of the type pp ^ cf) + X should be 



(1 + cos^ e) (2.66) 



8. There should be significant dependence on the beam and target polariza- 
tions in the reaction p-l-p p+p+(f), specifically, the cross section should 
be maximal when their polarizations are parallel. This effect should also 
be seen in the reaction p-l-d — > '^He-|-0, where there are recent indications 
of a substantial enhancement of 4> production relative to w production 
and the naive OZI rule. 

9. The constituent counting rules 

da{A + B^C + D) 
di 

suggest different behaviours for the reactions 7r + p-^(/),a; + nat large 
momentum transfers. Specifically, the cross section for n + p (j) + n 
should behave as s~^^ if it is dominated by production from the \uudss) 
Fock state, or as if it is dominated by gluonic production. In either 
case, the cfi/uj ratio should decrease at higher energies. 

10. The different production mechanisms for 07r and wtt, and for (jjinr and 
uj TTTT final states, could lead to different angular distributions, as may be 
seen experimentally. 

11. There should be a large enhancement in the Pontecorvo reaction 
p + d (f> + n. 



f{OcM) 



"CM 
fixnd 



nA+riB+nc+riD-'i 



(2.67) 



2. 6 Extension of the Model to A Production 
We finish this lecture by mentioning a couple of tests of these ideas 



122,123 



using 

data on A production in pp annihilation and elsewhere. 

The total cross section for pp — > AA and its angular distribution have 
been measured at LEAR by the PS 185 coUaborationlHfl. These have been 
described in terms of phenomenological models based on quarks and gluons, 
and on meson exchanges. The PS 185 dataE£fl indicate that the AA pairs are 
produced mainly in a spin-triplet state, as would be expected in our model 



(2.65), or in a gluon-exchange modelE^j, which produces a spin-triplet ss pair. 
Spin-triplet dominance can also be accommodated in a meson-exchange model 
123, with the appropriate combination of K and K* exchanges. Measurements 
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of the A depolarization Z3„„, i.e., the transfer to the final-state A of polarizaiipn 
from a polarized p target, may help distinguish between different modelsiif — 
Models which explain the EMC spin effect in terms of polarized gluonsEl"tJ 
would naturally expect that the gluon polarization is the same as the target p, 
and hence also that of the ss pair, and thus that of the A, so that Dnn > 0. On 
tlip other hand, the meson-exchange model has been used to predict Dnn < 
On the other hand, in our polarized-strangeness model we would expect 
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that the s in the prolun wave function would have negative polarization, and 



hence that Dnn < OliiJ. Thus the polarized-gluon and polarized-strangeness 



explanations of the EMC spin effect seem to predict opposite signs. 

This idea can be extended to A production in the target fragmentation 
region of dee p- ine lastic lepton-nucleon scattering. Briefly stated, the proposal 
made in Ref. [123| is that a polarized lepton beam (a i^, D, or a polarized fj, or e) 
couples to the nucleon target preferentially through a particular boson (W or 
7) polarization state, which then picks out preferentially a particular quark (or 
antiquark) polarization state. The target nucleon remnant "remembers" the 
spin that was removed, without the need for a polarized target. For example, 
if a positively-polarized u quark is removed from the proton wave function, 
it will tend to leave behind a negatively-polarized ss pair, which may lead-t^ 
negative polarization for A's produced in the target fragmentation regiont 
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as also seen in Fig. 2.12 





Figure 2.12: Diagrams which suggest that A baryons observed in the target fragmentation 
regions of deep-inelastic (a) P and (b) polarized /i+ collisions should be polarized. The solid 
triangles represent the spin states of the various particles. 



Just such an effect has been seen in the reaction ly + N^^^ + A + X 
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by the WA59 collaboration E£3. In the kinematic region 0.2 < xbj < 1, z < 0, 
they measure 

Pa = -0.85 ±0.19 (2.68) 



to be compared with the postdiction 

Pa = -0.94D 



(2.69) 



in our model, where D is an uncalculable dilution factor. The comparison 
between equations ( 2.67 ) and ( 2.68 ) above indicates that D = 0.9±0.2. Similar 
predictions can be made, for example, for A's produced in the NOMAD u 
experiment at CERN. 

The fact that, according to the WA59 experiment, there does not seem to 
be dilution by a large factor encourages the-extrapolation of these ideas to po- 
larized /i(e)"'' scattering, where we predictllH3 that in the target fragmentation 
region 

Pa = 0.7P^(e)i? (2.70) 

This should be observable in thejjxpposed HMC experiment at CERN, in 
the E665 experiment at FermilabllHS, and perhaps also in the HERMES O 
experiment at DESY. 
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3 Cosmological Spin-offs 



In this lecture we discuss, as an example of the relevance of polarized structure 
function measurements and the spin decomposition of the nucleon to other 
areas of physics, their applications to experimental searches for dark matter 
particles. This discussion is prefaced by a brief review of the motivations for 
such particles. 



3.1 How Much Dark Matter? 



Naturalness and inflation E£3 suggest that the density averaged over the uni- 
verse as a whole should be very close to the critical density, which marks the 
boundary between a universe that expands forever and one which eventually 
collapses, i.e. VL = p/pc —1- On the other hand, tlie-|matter we can see shining 



in stars, in dust, etc. amounts only to ~ 0.01 1123, as seen in Fig. |3.1| . The 
agreement between big bang nucleosynthesis calculations 11^ and the observed 
abundances of light elements suggests that ^baryons ^ 0.1, as also seen in 
Fig. Bjl This is to be compared with observations of rotation curves, which 
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suggest according to the virial theorem that the amount of matter in galactic 
haloes ilhaio — 0.1. A similar abundance of baryonic matters is suggested 
by observations of rich clusters of galaxies. 

Mathematically, the galactic haloes could in principle be purely baryqui 
although they seem unlikely to be made out of gas, dust or "snow balls" 
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However, there has recently been considerable interest in the possibility that 
haloes might be largely composed of "brown dwarfs" . There have been several 
searchesJbr such "failed stars" via microlensing of stars in the Large Magellanic 



Cloud which indicate that 

/ = 0.20t°:?i (3.1) 

be composed of brown dwarfs. This suggests that most of the dark matter 
present locally in our galactic halo 

Phaio = 0.3 GeV/cm' x 1.5°^! (3.2) 

is not baryonic in nature. Moreover, models of galaxy formation suggest that 



it is unlikely to be composed of massive neutrinosE^j. This leaves us only with 
the alternative of cold dark matter, whose detection involves in an essential 
way the spin decomposition of the nucleon, as we shall see later in this lecture. 

Before addressing in more detail the nature of the dark matter, we first 
comment on the age and Hubble expansion rate of the Universe, which have 
recently generated some controversy. Globular clusters seem to be at least 
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0.001 0.01 0.1 1 10 

Q. (density/critical) 



showing that there is no serious discrepancy 

10 - 

obtained from 



Figure 3.1: The (Q, Hq) plane, adapted from 
between the average measured value of Ho , f2 
This plot also shows the estimates of the present baryon densi 
visible features in the Universe, from Big Bang Nucleosynthesis 
All the indications are that i^taryons is 0.1, so that at least 90% of the matter in the Universe 

is non-baryonic dark matter. 



^baryons 

and from rich clusters. 



14 ± 3 Gyr old, and nucleocosmochronology suggests an age of 13 ± 3 
The question is whether these ages are compatible with current estimates of 
the Hubble-constant Hq km/sec/Mpc. Recent determinations of Hq may be 
combin ediH a to yield an estimate of 66± 13, which is shown on the vertical axis 
of Fig. 34, We see from this that there is no incompatibility between the .age 
of the Universe being above 10^° yr old and SI = 1 as wanted by inflationll2j. 
However, if this is indeed the case, at least 90% of the matter in the universe 
must be unseen dark matter. 
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3.2 Hot or Cold Dark Matter? 



These terms refer to whether the dark matter was relativistic or non-relativistic 
at the cosmological epoch when structures such as galaxies and clusters be- 
gan to form. Whether you favour hot or cold dark matter depends on your 
favourite theory of structure formation. If you believe that its origins lie in an 
approximately scale-invariant Gaussian random field of density perturbations, 
as suggested by inflationary modelsiHa, then you should favour cold dark mat- 
ter. This is because it enables perturbations to grow on all distance scales, 
whereas relativistic hot dark matter escapes from small-scaLe-^erturbations, 



whose growth via gravitational instabilities is thereby stunted t^. Thus galax- 
ies form later in a scenario based on Gaussian fluctuations and hot dark matter 
then they would in a scenario with cold dark matter. For this reason, the latter 
has commonly been regarded as the "standard model" of structure formation. 
Howeveu^if you believe that structures originated from seeds such as cosmic 



strings then you should prefer hot dark matter, because cold dark matter 



would then give too much power in perturbations on small distance scales. 




Figure 3.2: A compilation of data on the primordial perturbation spectrumt£il, compared 
with a cold dark matter simulation assuming an initially scale-invariant spectrum of Gaussian 

fluctuations. 
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Fig. 3.2 shows a compilationt^J of data oa-tlie power spectrum of astrophysical 



perturbations, as obtained from COBEEIHI and other observations of the cos 



mic microwave background radiation, and direct astronomical observations of 
galaxies and clusters. The solid line which does not quite pass through all the 
points is one calculated in the above-mentioned standard model of Gaussian 
fluctuations and cold dark matter. The discrepancies from this curve indicate 
that there is less perturbation power at small distances than would be expected 
in this theory, given the COBE normalisation at large distance scales. 

This and other observations have suggested that it may be necessary to 
modify the pure cold dark matter model. Several suggestions have been offered, 
including a non-zero cosmological constant and a deviation of the spectrum of 
Gaussian perturbations from scale invariance. However, the preferred scenario 
seems to be an admixture of hot dark matter together with the cold, resulting 
in the following cocktail recipe for the Universelilil: 

^cold — O-T, ri/iot — 0.2, Qharyons^O-^ (3-3) 

The way in which this scenario works is illustrated in Fig. [3.3|. Hot dark matter 




Scale 



Figure 3.3: Illustration how a 
perturbations seen by COBE 



««xed dark matter scenario Ela may reconcile the large-scale 
^'^^■^ with the relatively small magnitude of the perturbations 
seen at small scales. 



alone would give a spectrum of perturbations that dies out at small scales. 
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whereas hot dark matter does not. Combining the two, one can reconcile the 
relatively high COBE normalisation at large scales with the relatively small 
perturbations seen at small scales. 

The only realistic candidate that particle physicista^are able to offer for 



the hot dark matter is a neutrino weighing about 10 eVlil3, but there are many 
candidates for the cold dark matter, of which we discuss two in the rest of this 
lecture. 

3.3 Lightest Supersymmetric Particle 

This is our favourite candidate for cold dark matter Ei, and one for which 
the polarized structure function measurements are particularly relevant, as 
we shall see in the next section. The lightest supersymmetric particle (LSP) 
is expected to be stable in many models, and hence present in the Universe 
as a cosmological relic from the Big Bang. This is because supersymmetric 
particles possess a multiplicatively-conserved quantum number called R parity 
Efd, which takes the values +1 for all conventional particles and —1 for all 
their supersymmetric partners. Its conservation is a consequence of baryon 
and lepton number cancelation, since 

i?=(-l)3^+^+2S (3.4) 

There are three important consequences of R conservation: 

1. Sparticles should always be produced in pairs. 

2. Heavier sparticles should decay into lighter ones. 

3. The LSP should be stable, since it has no legal decay mode. 

If the LSP had electric charge or strong interactions, it would have con- 
densed into galaxies, stars and planets such as ours, where it could in principle 
be detected in searches for anomalous heavy isotopes. Electromagnetic or 
strong interactions would bind such an LSP deeply inside some conventional 
nuclear species. Various searches have excluded such bound LSPs above the 
abundances Eds 

4^^^^ ^ 10-1^ to 10-30 (3.5) 
n(protons) 

for 1 GeV ^ mj^sp ^ 10 TeV, which are far below the abundances calculated 
for such relic particles, which usually lie in the range 

rt (relics) 



n (protons) 



> 10-*" (em) to 10"'^ (strong) (3.6) 
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We conclude tiS that any supersymmetric relic LSP should be electromagnet- 
ically neutral and possess only weak interactions. Scandidates in the future 
sparticle data book include the sneutrino v of spin 0, some form of "neutralino" 
of spin 1/2, or the gravitino G of spin 3/2. The sneutrino is essentially excluded 
by the LEP experiments which measured the decay of the into invisible par- 
tirjpH which have counted the number of light neutrino species: 2.991 ± 0.0016 
which does not leave space for any sneutrino species weighing less than 
iM^, and by underground experiments to be discussed in the next section, 
which exclude a large range of heavier sneutrino masses. Since the gravitino 
is probably impossible to discover, and is aiwway theoretically disfavoured as 



143 



the LSP, we concentrate on the neutralino E 

The neutralino % is a mixture of the photino 7, the two neutral higgsi- 
nos Il\ 2 expected in the minimal supersymmetric extension of the Standard 
Model, and the zino Z. This is characterized essentially by three parameters, 
the unmixed gaugino T^i\|2^ the Higgs mixing parameter and the ratio of 
Higgs vacuum expectation values tan /?. The phenomenology of the lightest 
neutralino is quite complicated in general,pbut simplifies in the limit mi/2 0, 
where x is approximately a photino state 1113, and in the limit /i — > 0, where it 



is approximately a higgsino state. As seen in Fig. 3.4, experimental constraints 
fmm LEP and the Fermilab collider in fact exclude these two extreme limits 
11a, so that 

^ (10 to 20)GeV (3.7) 



Fig. 3.4 also indicates that there are generic domains of paramfiter space where 



the LSP may have an "interesting" cosmological relic densityEliI, namely 



0.1 < ?L^HI < 1 



(3.8) 



for some suitable choice of supersymmetric model parameters. Fig. 3.^iijsplays 
the calculated LSP density in a sampling of phenomenological modelsE£3 where 



we see that an interesting cosmological density is quite plausible for LSP masses 



20 GeV <, 



<, 300 GeV 



(3.9) 



The next problem is how to detect the LSP, and this is where polarized struc- 
ture function measurements may have a role to play. 



3.4 Searches for Neutralinos 

Several strategies to search for cosmological relic neutralinos have been pro- 
po,sed. One is to look for the products of their annihilations in our galactic halo 
Here the idea is that two self-conjugate x particles may find each other 
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Figure 3.4: The cosmological relic density of neutralinos xli^2 
in the range of interest to astrophysicists and cosmologists 



-»fl,y well (shaded regions) lie 
ila, namely 0.1 ^ ^ 1. 



while circulating in the halo, and have a one-night stand and annihilate each 



other: 



XX 



i£,qq, leading to a flux of stable particles such as p, e+,7, in 



the cosmic rays. Several experiments have searched for cosmic-ray antiprotons 
with the results shown in Fig. 3.5. At low energies there are only up- 



per limits, but there are several positive detections at higher energies, which 
are comparable withJJic flux expected from secondary production by primary 
matter cosmic raysE3. As also seen in Fig. 3.5, relic LSP annihilation in our 
galactic halo might produce an observable flux of low-energy cosmiii ray an- 
tiprotons somewhat below the present experimental upper limits t£fl. These 
may be interpreted as suggesting that 

Phalo 

(3.10) 



NA SA and the DOE have recently approved a satellite experiment called AMS 
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which should be able to improve signiflcantly the present upper limits 
on low-energy antiprotons, and may be able to start constraining significantly 
supersymmetric models. It is also possible to derive limits on such models fpaip 
the present experimental measurements of the cosmic-ray e+ and 7 fluxes 
but these are not yet very constraining. 



156 



53 




A second LSP detection strategy is to look for xx annihilation inside the 
Sun or Earth. Here the idea is that_a relic LSP wandering through the halo 
may pass through the Sun or EarthllHI, collide with some nucleus inside it, and 
thereby lose recoil energy. This could convert it from a hyperbolic orbit into an 
elliptic one, with a perihelion or perigee below the solar or terrestrial radius. 
If so, the initial capture would be followed by repeated scattering and energy 
loss, resulting in a quasi-isothermal distribution within the Sun or Earth. The 
resulting LSP population would grow indefinitely, a la Malthus, unless it were 
controlled either by emigration, namely evacuation from the surface, or by civil 
war, namely annihilation within the Sun oij-Earth. Evaporation is negligible for 



X particles weighing more than a few GeVt^, so the only hope is annihilation 
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The neutrinos produced by any such annihilation events would escape from the 
core, leading to a high-energy solar neutrino flux {E^, ^ 1 GeV). This could be 
detected either directly in an underground experiment, or indirectly via a flux 
of upward-going muons produced by neutrino collisions in the rock. 

The polarized structure function measurements enter in the estimate of 
the X capture rate, which enters in the following general formula E£3 for the 
neutrino flux: 



i?^=2.7xlO"V("^x/™p) 



cr{xP XP) 



Px 



/300 kms"^ 



10-4Ucm^ y V0.3GeVcm--/ \ Vy 

(3.11) 

where we have simplified to the case of capture by the Sun, where proton tar- 
gets predominate. Here / is a kinematic function, and are the local 
density and mean velocity of the halo LSPs, and represents factors asso- 
ciated with the neutrino interaction rate in the apparatus. The factor which 
interests us here is the elastic LSP-proton scattering cross section cr(xP ~> XP)- 
To see how the polarized struGtuxe functions enter into the estimation of 
the elastic scattering cross section E£j, first note that the LSP interacts with 
hadrons via an effective four-fermion interaction of the general form xXQQj 
which is mediated by the exchanges of m assive particles such as the Higgs 
bosons and squarks q, as seen in Fig. 3^. This four-fermion interaction is simi- 





q q 



Figure 3.6: The exchanges of massive particle such as the Z^, Higgs boson and squarks give 
rise to an effective four-fermion interaction between the neutralinos x ^'id quarks inside a 

proton target. 



lar in many ways to the original Fermi four-fermion weak interaction mediated 
by exchange. The matrix elements of the latter interaction between nu- 
cleon states are governed by the familiar /3-decay constant qa, which may be 
written in the form 

\gA\^Au-Ad (3.12) 

Analogous expressions in terms of the Aq exist for the spin-dependent part of 
the X"iiucleon scattering matrix element. For example, if the x were to be a 
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pure photino state, we would have 

4 11 
a„ = -Au+ -Aff+ -As (3.13) 
9 9 9 ^ ' 

It is amusing to note that this is exactly the same combinaiipn of the Aq that 



appears in charged-lepton scattering off a proton target In a sense, the 
EMC and its successors have been measuring 7-nucleon scattering! 

The EMC and subsequent measurements indicate values of the Aq that are 
rather different from those predicted by the nai've quark model, which means 
that the dark matter scattering matrix elements are also rather different, and 
hence also the upper limits on the halo dark matter density that can be deduced 
from a given search for high-energy neutrinos from the core of the Sun or the 
Earth. By now, as discussed in the first lecture, the determination of the Aq 
is quite precise, and, as discussed earlier in this lecture, we no longer believe 
that the LSP x can be a pure 7. Consider, for example, the plausible case 
of an essentially pure U{1) gauginoiSP B: its scattering matrix elements on 
protons and neutrons are given byEj 

17 ^ 5 
36^^+36^ 

(3.14) 

a„ ~ (Au <^ Ac?) 

We see that, in this case, the uncertainties from polarized structure mea- 
surements are likely to be much smaller than those from other components in 
equation ( ^.11| ). 

Some typical rate estimates for upward-going muons originating from high- 
energy solar neutrinos in a sampling of supersymmetric models are shown in 



.p - — Am + — (Ad + As) 



Fig. 3.7:: while-some models are already excluded by unsuccessful searches. 



most are not 
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We see in Fig. 3.8 that searches for solar signals usually 



constrain models more iJian searches for terrestrial signals, though this is not a 
model-independent factE£3. In the long run, it seems that a search for upward- 
going neutrino- induced JBiions with a 1 km^ detector could almost certainly 



detect LSP annihilationli£3, if most of the cold dark matter is indeed composed 
of LSPs. 

The third LSP seamh, strategy is to look directly for LSP scattering off 



nuclei in the laboratoryliiJ. It is easy to see that the typical recoil energy 



AE < m^v' 10 [j^) keV (3.15) 

deposited by elastic x-nucleus scattering would probably lie in the range of 
10 to 100 keV. The type of spin-dependent interaction, mediated by or q 
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Figure 3.7: The flux of upward-going muons expected from xxp*'pihilation inside the Sun 
in a samphng of supersymmetric models l£3- 



exchange, that we discussed in previous paragraphs is hkely to dominate for 
Ught nucleiE£3, whereas coherent spin-dependent interactions mediated by H 
and q exchange are hkely to dominate scattering off heavy nucleill£j. As we have 
abready discussed, the spin-dependent interactions on individual nucleons are 
controlled by the Ag: translating these into matrix elements for interactions 
on nuclei depends on the decomposition of the nuclear spin, which musi-Ja. 
studied using the shell model t£3 or some other theory of nuclear structurell^ 
The spin-independent interactions on individual nucleons are related to the 
different quark and gluon contributions to the nucleon mass, which is also an 
interesting phenomenological issue related to the 7r-nucleon cr-term discussed 
in Lecture II. Again, the issue of nuclear structure arises when one goes from 
the nucleon level to coherent scattering off a nuclear target. 

We will not discuss here the details of such nuclear calculations, but present 
in Figs, 



£ and [3.10 the results of a sampling of different supersymmetric 



models We see in Fig. 3.6 that the spin-independent contribution tends to 
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Figure 3.8: A comparison of the muon fluxes from the centre of t+teiSun and Earth, as found 
in a sampling of supersymmetric models t££l. 



dominate over the spin-dependent one in the case of Germanium, those th 
not universally true, and would not be the case for scattering off Fluorine 
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In Fig. [3.10 we plot the scattering rates off '^^Ge, where we see that there are 



many models in which more than 0.01 events/kg/day are expected, which may 
be observable. Thus the direct search for cold dark matter scattering in the 
laboratory may be a useful complement to the searches for supersymmetry at 
accelerators. 

3.5 Axions 



The axionE£3 is our second- favourite candidate for the cold dark matter. It 
was invented to guarantee conservation of P and CP in the strong interactions, 
which would otherwise be violated by the 9 parameter: 

JO-qcd 9 e'"'"" G^.Gp, (3.16) 
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Figure 3.9: A comparison of the spin-dependent and spin-independent interactioin I ' li tes of 
relic neutralinos x with Germanium in a sampling of supersymmetric models E£j- 



which is known experimentaUyB to be smaller than about 10~^. The 9 param- 
eter relaxes to zero in any extension of the Standard Model which contains the 
axion, a light pseudoscalar boson with mass and couplings to matter that are 
scaled inversely by the axion decay constant fa- 

Aqcd ruq TUf 1 

ma 7 > 9aff -J- : 5a77 J- (3-17) 

Ja J a Ja 

The fact that no axion has been seen in any accelerator experiment tells us 
that 

/a^lTeV (3.18) 

and hence that any axion must be associated with physics beyond the scale of 
the Standard Model. 
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Figure 3.10: Total scattering rates on Gernpmium in a sampling of supersymmetric 

models liia. 



Axions would have been produced in the early Universe in the form of 



slow-moving coherent waves that could const 
density of these waves has been estimated as 
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e cold dark matter. The relic 



0.6 X lO^'^cV 



7/6 



200McVy'^^^ 



A 



QCD 
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which is less than unity if 



fa 10^2 GeV 



(3.19) 



(3.20) 



In addition to |th^se coherent waves, there may also be axions radiated from 
cosmic strings 
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which would also be non-relativistic by now, and hence 
contribute to the relic axion density and strengthen the limit in equation (3.20). 

The fact that the Sun shines photons rather than axions, or, more ac- 
curately but less picturesquely, that the standard solar model describes most 
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data, implies the lower limit 
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/,;>10^GeV 



(3.21) 



This has been strengthened somewhat by unsuccessful searches for the axio- 
electric effect, in which an axion ionizes an atom. More stringent lower bounds 
on fa are provided by the agreementa^between theories of Red Giant and White 
Dwarf stars with the observations E£j : 



fa ;t 10^ GeV 



(3.22) 



Between equations (3.2C) and ( 3.22| ) there is an open window in which the axion 
could provide a relic density of interest to astrophysicists and cosmologists. 

Part of this window may be closed by the observations of the supernova 
SN1987a, which is where the polarized structure function measurements come 
into play. According to the standard theory of supernova collapse to form a 



neutron star, 99% of the binding energy released 
star escapes as neutrinos. This tiieory agrees 
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he collapse to the neutron 
with the_abservations of 



SN1987a made by the KamiokandeLLH and 1MB experiments tlrl, which means 



that most of the energy could not have been carried off by other invisible 
particles such as axions. 

Since the axion is a light pseudoscalar boson, its couplings to nuclear 
matter are related by a generalized Goldberger-Trciman relation to the cor- 
responding axial-current pi^trix elements, and these are in turn determined 



by the corresponding Aq 
individual nucleons that 
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Specifically, we find for the axion couplings to 



Cap = 2[-2.76 Au - 1.13 Ad + 0.89 As - cos 2/3 (Au - Ad - As)] 



Can = 2[-2.76Arf- 1.13Au + 0.89As-cos2/3(Ad- Au- As)] 



(3.23) 



Evaluating the Aq at a momentum seals around 1 GeV, as is appropriate in 
the core of a neutron star, we estimate E£I that 



Ca 



Ca 



(-3.9 ± 0.4) - (2.68 ± 0.06) cos 2/3 
(0.19 ±0.4) + (2.35 ±0.06) cos 2/3 



(3.24) 



which are plotted in Fig. 3.11 
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Figure 3.11: Axion couplings to the nucleon, as determined in Ref. jGst using polarized 

structure function data. 



As in the case of LSP scattering, the uncertainties associated with polar- 
ized structure function measurements are by now considerably smaller than 
other uncertainties, in this case particularly those associated with the nuclear 
equation of state. The total axion emission rate from the core of a neutron 
star is approximately proportional to the combination 



Ca 



0.8 {Ca 



0.5 C, 



ap 



(3.25) 



which is plotted in Fig. 3.12 , together with the associaled with the errors in 
the Aq. Estimating axion emission rates in this wayEIj, it ^eeps that part 

is underway 
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of the previous axion window is still open, and an experiment 
which should be able to detect halo axions if they occupy this window. 



3. 6 The Role of Spin in the Forthcoming Revolutions in Particle Physics and 
Cosmology 

Measurements of the spin decomposition may not only overturn our under- 
standing of nucleon structure, but may also help usher in the revolutions in 
particle physics and cosmology for which we all yearn. A "Standard Model" for 
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Figure 3.12: The total axion emission rate from the eore of a neutron star, proportional to 
the combination Can^ + 0.8 (Can + Cap)^ + 0.5 Cap^, as determined in Ref. fefj]. 



the formation of structure in the Universe is now emerging, in which Gaussian 
perturbations laid down during an inflationary epoch are ampUfied by gravita- 
tional instabilities fed by cold dark matter particles. The COBEjat|servations 
of primordial fluctuations in the microwave background spectrum^ 
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were per- 
haps the first observational evidence for this picture, just as the discovery of 
neutral currents heralded the establishment of the Standard Model of elemen- 
tary particles. Confirmation of this emerging picture of structure formation 
would be provided by the observation of one or more species of dark matter 
particle, such as massive neutrinos and LSPs or axions. The discovery of any 
of these would also cause a revolution in particle physics by taking us finally 
beyond the Standard Model, as well as telling us how galaxies and clusters 
were formed. As we have seen in this last lecture, polarized structure function 
measurements control the couplings to matter of both LSPs and axions, and 
could well play an important role in their detection. 



Long live the Spin Revolution 
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